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Abstract 

A U(2) flavor symmetry can successfully describe the charged fermion masses 
and mixings, and supress SUSY FCNC processes, making it a viable candidate for 
a theory of flavor. We show that a direct application of this U(2) flavor symmetry 
automatically predicts a mixing of 45° for z/„ => u s , where v s is a light, right-handed 
state. The introduction of an additional flavor symmetry acting on the right-handed 
neutrinos makes the model phenomenologically viable, explaining the solar neutrino 
deficit as well as the atmospheric neutrino anomaly, while giving a potential hot dark 
matter candidate and retaining the theory's predictivity in the quark sector. 
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1 Introduction 



The pattern and origin of the quark and lepton masses and mixings remains a challenging 
question for particle physics. Although a detailed description of this pattern requires a 
theory of flavor with a certain level of complexity, the gross features may be described 
simply in terms of a flavor symmetry and its sequential breaking. 

One simple flavor structure is motivated by four facts about flavor: 

• The quarks and leptons fall into three generations, ^1,2,3, each of which may even- 
tually have a unified description. 

• The top quark is sufficiently heavy, that any flavor symmetry which acts on it non- 
trivially must be strongly broken. 

• The masses of the two light generations imply a phenomenological description in 
terms of small dimensionless parameters, {e}. 

• In supersymmetric theories, flavor- changing and CP violating phenomena suggest 
that the squarks and sleptons of the first two generations are highly degnerate. 

It is attractive to infer that, at least at a phenomenological level, there is a non-Abelian 
flavor symmetry which divides the three generations according to 



The four facts listed above follow immediately from such a structure, with {e} identified 
as the small symmetry breaking parameters of the non-Abelian group. These control 
both the small values for quark masses and mixing angles, and also the small fractional 
non-degeneracies of the scalars of the first two generations. 

The Super-Kamiokande collaboration has provided strong evidence for an anomaly in 
the flux of atmospheric neutrinos, which may be interpreted as large angle oscillations of 
Ufj, predominantly either to v T or to u s , a singlet neutrino |T|. This observation provides a 
challenge to the non-Abelian 2 © 1 structure: 

• v T is expected to have a very different mass from that of z/ ejAt , and to only weakly 
mix with them. 

• If the atmospheric oscillation is — > u s , what is the identity of this new singlet 
state, why is it light, and how could it fit into the 2 © 1 structure? 

There are a variety of possible reactions to this challenge. One possibility is to drop 
the 2©1 idea; perhaps the CP and flavor violating problems of supersymmetry are solved 
by other means, or perhaps supersymmetry is not relevant to the weak scale. Another 
option is to retain the 2 © 1 structure for quarks, but not for leptons, where the flavor 
changing constraints are much weaker. 



2©1 : 



■ipa © Va, a = 1, 2. 
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In this paper we study theories based on the flavor group U(2), which immediately 
yields the structure (]l|), giving the 2 © 1 structure to both quarks and leptons 0. The 
masses and mixings of the charged fermions and scalars resulting fron U(2) have been 
studied in detail, and significant successes have been identified ||. We add a right- 
handed neutrino to each generation, and find that the symmetry structure of the neutrino 
mass matrix automatically chooses to be a pseudo-Dirac state coupled to one of the 
right-handed neutrinos, resulting in — > v s with a mixing angle close to 45°. 



2 U(2) Theories of Quark and Charged Lepton Masses. 

The most general U(2) effective Lagrangian for charged fermion masses, at leading order 
in the U(2) breaking fields, is 

C = fofoh + ( ^4> a ^ah + MS ab + A ab )ij b h) (2) 

where <p a is a doublet, S ab a symmetric triplet, A ab an antisymmetric singlet of U(2), and 
h are Higgs doublets. Coupling constants have been omitted, and M is a flavor physics 
mass scale. An entire generation is represented by ip, so that each operator contains terms 
in up, down and charged lepton sectors, but unification is not assumed. For example, this 
theory follows from a renormalizable Froggatt-Nielsen model on integrating out a single 
heavy vector U(2) doublet of mass M (see the second of [[J). 

The hierarchical pattern of masses and mixings for charged fermions is generated by 
breaking U(2) first to U(l) with vevs <f) 2 , S 22 ~ eM, and then breaking U(l) via the vev 
A 12 fa e'M. The symmetry breaking 

U(2) =4 U(l) 4 1 (3) 
produces the Yukawa coupling textures 




M LI , r \ c> e c | (J) 



3 General Effective Theory of Neutrino Masses. 

Without right-handed neutrinos, the most general U (2) effective Lagrangian for neutrino 
masses, linear in £7(2) breaking fields, is 

C e ff = Tjhhhh + -L (l^ a l a hh + l a S ab l b hh) . (5) 
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where l a , I3 are lepton doublets. The term l a A ab lbhh vanishes by symmetry; hence the 
above vevs give the neutrino mass texture 



Mll = 77 e e . (6) 




so that the lightest neutrino is masslessj] The mixing angle for — > v T oscillations, 6^ T , 
is of order e — the same order as mixing of the quarks of the two heavier generations, 
- and is much too small to explain the atmospheric neutrino fluxes. However, in theories 
with flavor symmetries, the seesaw mechanism typically does not yield the most general 
neutrino mass matrix in the low energy effective theory. This apparent problem requires 
that we look more closely at the full theory, including the right-handed neutrinos. 



4 The Seesaw Mechanism: A Single Light vr 

Adding three right-handed neutrinos to the theory, N a + Ns, the texture for the Majorana 
mass matrix is: 

/0 0\ 

M RR = M e e . (7) 
\0 e l) 

with the 12 and 21 entries again vanishing by symmetry. In supersymmetric theories the 
zero eigenvalue is not lifted at higher order in the flavor symmetry breaking. This presents 
a problem for the 3x3 seesaw mechanism in U(2) theories, since Mll = MlrM rr M-[ r 
and M RR cannot be inverted. 

One approach || is to allow further flavor symmetry breaking vevs, for example (f) 1 7^ 0, 
so that Mrr has no zero eigenvalues. Remarkably, taking 1 /M pa e', the seesaw gives 
6^ pa 1, as needed for the atmospheric neutrino anomaly. On the other hand, this 
pattern of neutrino masses cannot explain the solar neutrino fluxes, and the additional 
flavor breaking vevs remove two of the highly successful mass relation predictions of the 
quark sector. 

In this paper we keep the minimal U(2) symmetry breaking vevs and pursue the 
consequences of the light N e state which results from (0). The singular nature of M RR is 
not a problem; it is an indication that iV e cannot be integrated out of the theory. However, 
Af- and do acquire large masses, and when they are integrated out of the theory the 
low energy 4x4 neutrino mass matrix is: 

( \ 

M LL e'v 



\0 e'v J 

1 Including operators higher order in the U(2) breaking fields, the lightest neutrino remains massless in 
a supersymmetric theory, but not in the non-supersymmetric case , where operators such as l a A ab '(fnlshh 



M (4) 



(8) 



occur. 
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where Mll is a 3 x 3 matrix in the (v a , v 3 ) space, determined from seesawing out the two 
heavy right-handed states, and has one zero eigenvalue. 

Because the N e — u^ mixing is weak scale, while all other couplings to are suppressed, 
N e and are maximally mixed. Thus, we note that a direct application of the U(2) theory 
to the neutrino sector predicts a 45° mixing between and vj 

There is a significant phenomenological difficulty with this model. The mass of the 
N e — pseudo-Dirac state is of order e'v. Using a value for e' extracted from an analysis 
of the charged lepton sector, this is of order 1 GeV, well in excess of the 170 keV limit 
obtained from direct searches. One simple solution is to restrict the couplings of the 
right-handed neutrinos by an additional U(1)n approximate flavor symmetry. Each N 
field carries N charge +1, while the symmetry is broken by a field with charge -1, leading 
to a small dimensionless breaking parameter ejv- The entries in the neutrino mass matrices 
receive further suppressions 

M LR ^e N M LR M RR ^t 2 N M RR (9) 

which, for the 4x4 light neutrino matrix, simply leads to the replacement e'v =>- e^e'v 
in the N e — entry, giving 



M< 4 » = 



I e M M M u i 

/ v 2 v 2 v 2 / 

€ M € M e M e ^ ev 

I v 2 v 2 v 2 n 

fc M t M M U 

V e N e'v J 



(10) 



It is understood that all entries have unknown 0(1) coefficients. 

Note that Mll is unchanged. There is a simple reason for this. If we modify our 
right-handed couplings by the replacements Mlr — > MlrT, M rr — > T T M RR T, where T 
is any diagonal matrix, then 

M LL => Ml R T{T t M rr T)-\Ml R T) t = M LL - (H) 

It is interesting that the observed value of 8m 2 Q can give the appearance that right-handed 
neutrinos receive GUT-scale masses, while their masses are in fact much lower. 

If the iV e — entry dominates the mass of i.e. if ^> jfe, this 4x4 matrix 
splits approximately into two 2x2 matrices, and maximal mixing is preserved. One 2x2 
matrix describes the pseudo-Dirac state 



( e i e Nt'v 
\e N e'v 



(12) 



while v e =3- u T mixing is described by 




(13) 
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(1) Heavy states 

(2) Light states 


ve N e' - ve N e' + 

e M M 


£e N ee> 45° 



Table V. General Theory: the masses, mixings, and splittings of the two sets of neutrinos. 

The resulting masses and mixings are given in Table [I]. 

Since e and e' are determined by the charged fermion masses, in the neutrino sector 
there are two free parameters, and M, which describe five important observables: 9®, 
atm , 5m 2 & , 5m 2 atm and m u , the mass of the pseudo-Dirac muon neutrino. However, the 
various predictions of the theory have varying levels of certainty. Because there are a large 
number of order one constants in the original formulation of the theory, we can end up 
with a prediction which has a coefficient of a product of some number of these quantities. 
To assess the level of certainty, we will include a quantity i, which we term the "stability 
index" of the prediction, which is simply the power of unknown order one coefficients 
appearing in the prediction. 

Two of the three resulting predictions are the mixing angles 

sin0 «e' [z = 4], 9 atm = 45° [i = 0] . (14) 

The postdiction of a maximal mixing angle for atmospheric oscillations is an important 
consequence of the U(2) theory. The value of e' extracted from the charged fermion sector 
is 0.004, within an order of magnitude of the central value 9® = 0.037 of the recent BP98 
fit to the solar data, and within a factor of 4 of the minimal acceptable value of 0.016 
||. Such a discrepancy is not a great concern, as we gain a comparable contribution from 
the charged lepton matrix. Furthermore, the prediction of 9 Q involves the fourth power 
of unknown order one coefficients, thus % = 4, and is somewhat uncertain. 

The relevant mass splitting for the u e — > v T oscillations occuring in the sun is 

SK«(ff. (15) 

While this is not a prediction of the theory, it is intriguing, as has been noticed elsewhere 
in other contexts, that if M is taken close to the scale of coupling constant unification, 
StUq ~ 10~ 5 eV 2 , in the right range for either small or large angle MSW oscillations. 

The final free parameter is fixed by the observed mass splitting for atmospheric 
oscillations 

i 

tmltm ~ ee ' eN M ~ ee ' eNV V 5m ® ( 16 ) 
giving ~ 10~ 8 — the U{1)n symmetry is broken only very weakly. 
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The final prediction is for the mass of the heavy pseudo-Dirac u^N e state: 

A 777 ^ 

m v « e'e N v « p!!L « 10°' 4 eV - 10 2 eV, [i = 4] (17) 
ey]8m% 

where the given spread in mass is due to uncertainty in 5m 2 atm and StUq. While it is 
tempting to interpret this as a good candidate for hot dark matter, we will see later that 
KARMEN places stringent limits on the acceptable values of m v . 



5 A Variant Theory 

A variation on this breaking structure was explored in a particular model (see the second 
of |0]), and it is interesting to explore whether this same approach for neutrino masses 
can work within that model. In this variation, there is no S ab field present, and the RR 
and LR masses are given by 



M, 



LR 






e' 






-e' 







M RR 





e 


; 






(18) 



generating a light 4x4 mass matrix 



/ iLvL 

1 e 2 M 
c M 
e M 





/ v 2 
fc M 



e— 

e N e'v 



' M 

t M 

M 





e N e'v 





(19) 



This matrix is problematic, because the 2x2 submatrix for the atmospheric neutrinos 
does not contain a splitting term. Of course, a splitting would be generated through 
interactions with the other left-handed states, we estimate 



Consequently our atmospheric splitting is 



m. 



'AT 



(20) 



(21) 



5rriQ, this would predict 5m 2 atm < <5m|, which is unacceptable. One 
simple solution is to allow the appearance of the operators 



Since we have {jj) 2 



{—Y<t> a <t> b N a N b M GUT 



(22) 
(23) 
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"llight TTlfieavy 


Sm 2 9 mix 


(1) Heavy states 

(2) Light states 


n,a J V 2 J j_ 2 V 2 

ve N e - e m ve N e + e m 

fe'\2v 2 v 2 
\ e 1 M M 


v 3 c 2 1 A pro 

(v 2 \2 e' 
v M > t 



Table 2: Without S field: The masses, mixings, and splittings of the two sets of neutrinos. 



The inclusion of one or both of these operators in our Lagrangian has the same effect 

-(4) ^ 2vi 
MM ~ M 



on our final mass matrix, inducing Mf$ ^e 2 ^ and yielding the 2x2 sub matrix 



describing the pseudo-Dirac state, while v e =>■ v T mixing is now described by 

The resulting masses and mixings are given in table |2|. 
The mixing angles in this variation are predicted to be 



sin(0 )«^ [i = 5], 6 atm = 45° [i = 0] (26) 

As the pseudo-Dirac muon neutrino is still present, the atmospheric angle is unchanged. 
However, the solar angle is changed somewhat. We should note that values for e and 
e' extracted for a fit of this model are different than for those of the previous model. 
Using values from fits in the charged fermion sector, we have e ~ 0.03 and e' ~ 5 x 10~ 4 
or e' ~ 2.4 x 10 -4 (depending on certain signs), yielding 9 & « 0(1.5 x 10~ 2 ). Given 
the number of 0(1) parameters involved, this is again quite consistent with the BP98 
small-angle MSW solution. 

The solar splitting scale is unchanged, while the atmospheric splitting is further sur- 
pressed by a factor of e. 



$ m at m « t'^Nve \J5m% (27) 

We fit this splitting again with the free parameter ~ 10 -6 — 10~ 7 . The resulting muon 
neutrino mass is then 

()Tn 2 

mu ~ « 10 L7 eV - 10 35 eV [i = 5] (28) 

8m% 
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Thus, while the explanations of the solar and atmospheric neutrinos remain, the neu- 
trino becomes potentially dangerous in its cosmological implications. However, given the 
large stability index of this prediction, there are large uncertainties in the prediction for 
its mass. 



6 KARMEN and LSND 

The presence of an additional sterile state makes it possible that a signal would be seen in 
short baseline —>■ v e oscillations, such as has been reported at LSND M. An estimate of 
the LSND mixing angle from the neutrino sector gives ee' 5m 2 & / 5m 2 atm , a very small result. 
Hence, this mixing originates from the charged lepton sector 



Olsnd = J— [i = 0]. (29) 

V m M 

The precise predictions for 1 — 2 mixing angles in the charged sector is an essential feature 
of the U(2) flavor symmetry. In the quark sector it is highly successful. In the lepton 
sector, Qlsnd = is om y useful if the neutrino mixing is either predicted or small, 

as in this theory. Recently, the KARMEN experiment has placed limits on the allowed 
region for such oscillations, giving a limit m v < 0.6 eV ||. While the prediction for m v 
has a large stability index in both the general theory as well as the variant theory, because 
the initial range for m v is so high in the variant theory, it is disfavored by this bound. 

The general theory is much safer, however. As we discuss in the appendix, the uncer- 
tainty due to order one coefficients would allow it to satisfy the KARMEN bound. Such 
a result would likely coincide with higher values of SttLq and lower values of 8m 2 atm . 



7 Astrophysical and Cosmological Implications 

There are three important cosmological implications of our theory. 

1. We predict a small, but potentially significant amount of neutrino hot dark matter. 
The KARMEN bound limits us to a 0.6 eV neutrino, but because there are two massive 
states, it is still within the interesting region for HDM. 

2. We predict abundances for light nuclei resulting from four light neutrino species. While 
newer data suggest D/H ratios lie in the low end of the range previously thought, and 
thus N v < 4, this is still an open question. 

3. There may be two further singlet neutrino states, dominantly and N T , at or below 
the weak scale. Successful nucleosynthesis requires that they decay before the era of 
nucleosynthesis. Because the mass eigenstates are slightly left-handed, the primary decay 
mode will be through the process shown in figure |]. This is similar to muon decay, which 
we use as a benchmark. For the lighter of the two states, we estimate its lifetime to be 
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5m% , 



m, 



The mass of this particle is 



m N u 



1 {S™ 2 atm) 



(5m|) 3 / 2 



for the general theory and 



m N u 



1 (5m 



2 ^ 
atrrii 



(5m%fl 2 



[i = 12] 



[i « 11 - 16] 



(30) 



(31) 



(32) 



in the variant theory. The stability index is approximate because it involves sums of order 
one coefficients of different powers. Furthermore, % will change depending on which of (|23| ) 
are included. 

The more dangerous case, the general theory, then has a mass O(100MeV) and thus a 
lifetime ~ 10 3 s, which is far too long to be acceptable. However, because the lifetime 
has a fifth power dependence on the mass, and because the prediction for the mass has 
index 12, deviations in the order one quantities could very easily push the lifetime down 
to an acceptable level. As we explore in the appendix, even conservatively we can only 
reasonably estimate the mass of this particle to be in the range (17MeV, 40GeV), which 
means that the lifetime could easily be 10 _9 s, without even beginning to push the limits 
of the order one quantities. The details are presented in the appendix. 



8 Models 

The theory described in this paper has a low energy effective Lagrangian of @) for charged 
fermion masses, while the neutrino masses arise from the U(2) x U(1)n effective La- 
grangian 

W = jfN 3 l 3 h + ^ ( N 3 <j> a l a h + h<P a N a h + N a (S ab + A ab )l b h 

+ & (N 3 N 3 M + N 3 N a ct> a M + N a N b S ab M) (33) 

where A^ 3 and N a have U(1)n charges +1, while 4>n has U{1)n charged —1. The field 
(f>N gets a vev, breaking U(1)n and establishing an overall scale for these coefficients: 
< ^ > = e^y. This effective theory can result from a renormalizable model by integrating 
out heavy states, both singlet and doublet under U(2), in the Froggatt-Nielsen mechanism. 

This symmetry structure on the right-handed singlet sector is far from unique. Another 
possibility is for N a to carry U(1)n charge, while A^ 3 is neutral under U(1)n- This has no 
effect on any of our predictions, since the form of (|H]) for the light neutrino mass matrix 
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Figure 1: Principal decay mode for N^. 



10 



is unchanged. The only change is that A3 has a mass of the order of the unification scale 
M rather than of order e 2 N M . 

Another possible symmetry structure for the theory is U(2)$ x U(2) N , where U{2)^ 
acts as usual on all the matter with non-trivial SU(3) x SU(2) x U(l) quantum numbers, 
while U(2)n acts only on the three right-handed neutrinos, with A 3 a singlet and iV a a 
doublet. The matrix Mrr now has the form of (0), and arises from the renormalizable 
interactions 



with vevs for S 22 and <j) 2 being of order eM and breaking U{2)^ — > U(1)n. The interac- 
tions for M LR are 



where R aA transforms as a (2,2). The vev for R 22 is also of order eM, since this is the 
scale of breaking of U(2)^ x U{2) N — > U(l)^ x £7(1)^. The breaking scale for U{1)^, is 
e'M, so the vev of R 12 takes this value. On the other hand, U(1)n is broken by R 21 . We 
choose this scale to be smaller by a factor of ejv, < R 21 >~ e^e'M, giving 



Integrating out the heavy states A 2 and A3, which now have masses of order the unification 
scale, this theory now reproduces ( |T0D for the mass matrix of the four light neutrinos. 
The common features of these models, which are inherent to our scheme, are: 

• There is a U(2) symmetry, which acts on the known matter as ip3®ip a , and is broken 
sequentially at scale eM and e'M. 

• A U(2) symmetry also acts on the three right-handed neutrinos with A3 a singlet 
and A^ 12 a doublet. This U(2), together with the symmetry of the Majorana mass, 
implies that Ni does not have a Majorana mass and becomes a fourth light neutrino. 

• There is an addition to the flavor group, beyond the U(2) which acts on ip. At least 
part of this additional flavor symmetry is broken at a scale very much less than M, 
leading to a small Dirac mass coupling of v^,N e . Such a small symmetry breaking 
scale could be generated by the logarithmic evolution of a scalar m 2 term. 

9 Conclusions 

There are several theories with sterile neutrinos ||, [ID], [Tl[ some of which have 4x4 textures 
that split into two 2x2 matrices. Such theories provide a simple picture for atmospheric 



Wrr = MN 3 N 3 + N^ A N A + N A S AB N B 



(34) 



W LR = l 3 N 3 h + — {l a ej) a N z h + h<j) A N A h + l a R aA N A h) 



(35) 




(36) 
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oscillations via — > v s , and solar oscillations via u e — > z/ T , with « ^ 10~ 5 eV 2 for 
M M uni f. However, theories of this kind typically do not provide an understanding for 
several key points: 

• Why is the Majorana mass of the singlet state v s small, allowing v s in the low energy 
theory? 

• Why does v s mix with rather than with v e or v T 1 

• Why is the v s — state pseudo-Dirac, leading to 45° mixing? 

• How can this extended neutrino sector be combined with the pattern of charged 
quark and lepton masses in a complete theory of flavor? 

• What determines the large number of free parameters in the neutrino sector? 

In the theory presented here, all these questions are answered: the key tool is the 
£7(2) flavor symmetry, motivated several years ago by the charged fermion masses and 
the supersymmetric flavor problem. The simplest pattern of U(2) symmetry breaking 
consistent with the charged fermion masses does not allow a Majorana mass for one of 
the three right-handed neutrinos. Furthermore, it is precisely this right-handed state that 
has a Dirac coupling to but not to v e or v T , guaranteeing that is pseudo-Dirac with 
a 45° mixing angle. 

Our theory provides a unified description of both charged fermion and neutrino masses, 
in terms of just three small symmetry breaking parameters and a set of order unity coef- 
ficients. Some predictions, such as = \jm u lm c and 9 atm = 45° are independent 
of the order unity coefficients and are precise. Other predictions, such as jV^I ~ m s /m& 
and #0 ~ ^Jm e m^/m 2 involve the order unity coefficients and are approximate. In the 
appendix we have introduced the "stability index" which attempts to quantify the uncer- 
tainty in such predictions according to the power of the unknown order unity coefficients 
appearing in the prediction. 

There is one further free parameter of the theory — the overall mass scale M setting 
the normalization of the right handed Majorana mass matrix. If M is taken to be the 
scale of coupling constant unification <5m ^ 10~ 5 eV 2 . 

The value of Sm 2 tm is not predicted — this is the largest deficiency of the theory. It can 
be described by a very small flavor symmetry breaking parameter. Once this parameter 
is set by the observed value of 5m 2 atm) it can be used to predict the approximate mass 
range of the pseudo-Dirac to be in the range 10 a4 — 10 2 eV, with significant additional 
uncertainty due to order one coefficients. This, even with the KARMEN bound, allows 
for a neutrino of cosmological interest with J2i m v % ~ 1 e V. Such a neutrino could be 
seen at short baseline experiments, and may have already been seen by LSND. Searching 
for — > z/ e , with sin 2 (2$) = 2 x 10~ 2 , below the current limit of 5m 2 is an important 
experiment for the U(2) theory, since it is this prediction which differentiates U(2) from 
several other theories with a light singlet neutrino. 
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Experiment 

Present solar v exp. 

SNO 

Borexino 

KAMLAND 

LSND,KARMEN 

K2K 

MINOS, ICARUS 
Atmospheric v exp. 



Mode 



V T 
V T 

N e 
N e 
N P 



Signal 

All data consistent with 2-flavor MSW 
Confirm SK measurement of B8. Measure 4^ 7^ 1 
Consistent with small-angle 2-flavor MSW 
No signal 

sin 2 (20) = 2 x 1(T 2 

disappearance. No e appearance 

disappearance. No r appearance. 
Confirm 2 flavor z/„ — > v s with 45° mixing. 



Table 3: Experimental signals. 



Predictions of the theory for experiments sensitive to neutrino oscillations are listed 
in table |3|. We expect a small angle MSW solution to the solar neutrino anomaly, through 
a v e =^> v T oscillation. The atmospheric neutrino anomaly is from v s . This will be 

distinguishable from => v T through a number of means: LBL experiments will see 
disappearance, but no v e or v T appearance. Improved statistics from Super-Kamiokande 
will be useful in distinguishing v T and =>■ u s , for example via inclusive studies of 

multi-ring events |fP2|. 



A "Formalism" of the Stability Index 

It is difficult to establish a formalism for the stability index, because it involves an in- 
herently ill-defined quantity, namely, what constitutes an order one quantity. However, 
the potential instability of various predictions to variations in these order one parameters 
makes some attempt to quantify this necessary. Such a quantification should be relatively 
insensitive to what precisely constitutes an "order one quantity" . 
Therefore, we demand the following quantities of the index: 

• An "order one" quantity should be defined as a quantity x with some probability 
distribution P(x) to occur in an interval about 1. For reasons that will become clear later, 
it will be useful to consider instead the quantity P(y), where x = 10 y . 

• This distribution should be "sensible", namely 

1. P(x) should be an even function in Log(x), that is, P(y) is even in y. 

2. P(y) should achieve its maximum value at 0. 
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3. P(y) should have a spread characterized by its variance, 



■oo 



y 2 P(y) 



— oo 



the variance then quantifying what "order one" is numerically. 

4. A product of two sensible distributions, correlated or uncorrelated, should be 
sensible. 

• The index should have similar implications regardless of P(y), so long as it is sen- 



• The definition of P(y) should be the only necessary input. 

We shall explore the motivation for these assumptions and will shortly see that the 
presented index nearly meets the requirements, and with minor modifications can meet 
them enitrely. 

We assume that the expectation value of x, and of any products of x, is unity. It 
follows immediately that P(y) should be even in y. We do not have strong arguments in 
favor of this assumption, and if it were relaxed, the formalism could be suitably modified. 

For instance, consider the seemingly sensible distribution 



which has been normalized to give total probability 1. The expectation value of a product 
of n uncorrelated variables with such a distribution would be 



Such a numerical pile-up of the central value of a product of order unity coefficients is 
excluded by our assumption. 

What constitutes a "sensible" distribution is, of course, a judgement call. Examples 
of what we consider sensible distributions would be 

• Flat distributions taking on the value 1/a from —a/2 to a/2 

• Exponential distributions with standard deviation a 

• Linearly decreasing distributions of the form 



sible. 





(37) 





otherwise. 
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In fact, it can be shown that the last case it just the product of two uncorrelated quantities 
of the first type. 

In all of these cases, the next moment (x A ) is irrelevant in quantifying the likelihood 
of the variable being within a particular region about zero. Requirement 3 is then simply 
a statement that a sensible distribution should simply have one quantity, its variance, to 
determine how confident we are that the variable is within that region. This will then 
allow us to be more confident in deducing the significance of the variance of some product. 

This being stated, we can actually go about constructing some approximation of con- 
fidence intervals. The ability to describe the distribution of one variable by its variance is 
useful in allowing us to calculate the variances for higher products. We begin by writing 
the formal expression for the probability distribution of n uncorrelated variables Xj = lO Vi 
with probability distributions Y\ i P(y i ). We have 

P(z) = J JryiRPiiyiMz - J» (38) 

i i 

This expression is tedious to calculate for given P(y), particularly for large n. However, 
its variance is a relatively simply calculation. 

v\ = J dzP{z)z 2 = J dzd n y{Y,yi)\J{Pi{yi))^-Y,yi) ( 39 ) 

" i i i 

Expanding the squared term we find terms 

fdzdryy i y s (^P i (y i ))8(z-'£y i ) = [ ^ (40) 

giving 

v 2 =Yv 2 (41) 

^zuncorr u i V^^J 

i 

For n correlated variables, a similar calculations yields 



(42) 

where v 2 , is the variance of the original variable. 

Thus, a product of n correlated order one quantities is far more unstable than a product 
of n uncorrelated order one quantities. Simply counting the total number of order one 
coefficients is not sufficient. Thus we will refer to a product of the form 

IK* (43) 

i 

as having index (J2i n i) of type (ni,7i2, ...,n m ). If some of the rii are repeated, we use the 
shorthand of writing n\ if n is repeated j times. We assume all order one quantities have 
the same distribution. A product of type (ni,n 2 , ...,n m ), has variance 
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4h) = ffln? ( 44 ) 

This works extremely well for products of order one coefficients. However, a sum of 
order one coefficients is not necessarily order one. In these cases, it is usually best to 
perform a Monte Carlo to determine the distribution. 

B Sensible distributions 

To characterize the probability of a general product to be within a certain region about 1, 
it is necessary to explore the particular forms of various distributions. We consider three 
reasonable distributions to be i) the flat distribution, ii) the Gaussian distribution, and 
iii) the linearly decreasing distribution. 

A product of two equal width flat distributions yields a linear distribution, so we need 
only consider the flat and Gaussian cases. Gaussian distributions are well understood: 
products of variables with Gaussian P(y) functions are again Gaussian, allowing standard 
statistical techniques to be applied. 

Products of flat distributions very quickly become characterized by Gaussian distri- 
butions. We have performed explicit Monte Carlos for n — 1,2, 3, 5, 7, 9 uncorrelated 
variables. Even by n = 2 the Gaussian approximation is good, and for n > 3 it is very 
good. We thus believe it is reasonable to simply use Gaussian distributions, making a 
statistical interpretation of the variances simple. 

For a standard, we propose using a distribution with variance v = y^, which cor- 
responds to the variance of a flat distribution for — | < y < |. Changing the width of 
such a distribution from 1 to a would amount to multiplying this variance by a. Such 
generally mild sensitivity of the index to variations in the initial distribution is one its 
desirable qualities. We can then take "l-i>" and "2-t>" regions with \y\ < v and \y\ < 2v, 
respectively. As should be clear, these should not be interpreted as the precise 67% and 
95% 1 and 2-cr regions, because a is not precisely defined. They are simply regions of 
medium and strong confidence, respectively. 

As an example, consider a prediction with an unknown coefficient of order one quan- 
tities of the form x\x2X%. We say this has index 2 + 1 + 1 = 4 of type (2, 1, 1), which 
we will write in shorthand as (2, l 2 ). Assuming the standard variance given above, this 

coefficient has variance v = y / 22+ ^ +1 — = ^/|. Thus, we can have medium confidence that 
the prediction for x is known within a factor of 10 v = 5, and strong confidence the the 
prediction is within a factor of I0 2v = 25. 

We can also see that this reduces to the expected prediction in the case of a variable 
of index 1 of type (1). It will have variance v = y-^ which gives medium confidence that 
the prediction is known within a factor of 1.9, and strong confidence it is known within 
a factor of 3.8. This is a good consistency check that the index predicts what we would 
expect in the case of a single order one coefficient. 
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quantity 


type 


Sign 


Range of Med 


Range of Strong 




(*=approx type) 


convention 


Confidence 


Confidence 


a 





n/a 


? exact 

4 


exact 


e Q 


4* 


+ 


(0.002,0.03) 


(0.0006,0.1) 


e e 


4* 




(0.005,0.012) 


(0.0001,0.06) 


m v 


(2, 1 2 ) 


none 


(2, 50)eV 


(0.4,250)eV 


m NL 


(4,3,2 2 ,1) 


none 


(17MeV,40GeV) 


(300keV, 1.9TeV) 


m N[I 


(4, 2 4 ) 


none 


(470MeV, 870GeV) 


(HMeV, 37TeV) 



Table 4: General Theory: uncertainties in predicitions. The regions listed here are simply 
for the uncertainty due to order one coefficients. Additional error due to uncertainty in 
input quantities, in particular in m v , can also be significant. 

C Reassessing the uncertainties in the U{2) neutrino 
model 

In lieu of the preceding analysis, we address the index type of the predictions already 
presented, and thus assess strong and medium confidence regions of each prediciton. We 
list all uncertainties for the general theory in table |4]. 

In the general theory with the S-field, the atmospheric mixing angle is completely 
stable, while the solar angle is of approximate type (l 4 ) = (1, 1, 1, 1). However, it involves 
a sum of order one coefficients, motivating the use of Monte Carlos. Since a sum is 
involved, the relative sign of the order one quantities becomes relevant, and we list those 
cases seperately. These Monte Carlos allow us to claim that we have medium confidence 
that # lies within (0.002,0.03) and strong confidence that it lies within (0.0001,0.1), 
giving large overlap of the BP98 region. 

The mass of the pseudo-Dirac neutrino has stability index 5 of type (2, l 3 ), giving 
a medium confidence to know this within a factor of 5, and strong confidence within a 
factor of 25. Given the uncertainty in 5m 2 atm and SitLq, which determine the prediction, 
m u could conceivably be as low as O.leV. 

The masses of the right-handed states are not known so well. The mass prediction is, 
for the heavier state, of type (4, 2 4 ), and, for the lighter state, of type (4, 3,2 2 , 1). This 
would give medium confidence to know the masses at factors of 43 and 48, and strong 
confidence at factors of 1800 and 2300, respectively. The cosmological implications of 
these neutrinos are very uncertain, given that the lighter could be well over a TeV in 
mass. 
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Without the S field, certain uncertainties change. The precise nature of the changes 
depends on which splitting operators are included and what sign convention is taken. 
Because of the large number of permutations, we list only the basic results. The atmo- 
spheric angle is, as expected, completely certain. The solar angle becomes slightly more 
uncertain, but still overlaps BP98 well. The heaviest two righthanded masses typically 
become less certain by a factor of roughly 100, but the uncertainty is so large that the 
phenomenological predictions remain the same. The only dramatic difference in the vari- 
ant theory is that has a medium confidence region on its mass of (24eV, 4keV), and a 
strong confidence region of (2eV, 48keV). Including the uncertainties in the input quan- 
tities, the mass could be as low as OAeV, which escapes the KARMEN bound, although 
narrowly. 
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